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Abstract 

The quasiclassical correction to the Mohere's formula for multiple scattering is derived. The 
consideration is based on the scattering amplitude, obtained with the first quasiclassical correction 
taken into account for arbitrary localized but not spherically symmetric potential. Unlike the 
leading term, the correction to the Moliere's formula contains the target density n and thickness L 
not only in the combination uL (areal density). Therefore, this correction can be reffered to as the 
bulk density correction. It turns out that the bulk density correction is small even for high density. 
This result explains the wide region of applicability of the Moliere's formula. 
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I. INTRODUCTION 



Multiple scattering of high-energy particles in matter is a process which plays an impor- 
tant role in the experimental physics. The basis of the theoretical description of this process 
goes back to Refs. . Further development of the theory of multiple scattering was 

performed in numerous publications, see, e.g., Refs. a M and references therein. Detailed 

nn 

experimental investigation of multiple scattering has also been performed, see Refs.p, 18|. 

The celebrated Moliere's formula describes the angular distribution ^ for small-angle 
scattering. It was shown by Bethe in Ref. jj] that the most simple way to derive this 
formula is to use the transport equation. As a consequence of this equation, the angular 
distribution ^ depends on the thickness L and the density n only in the combination nL, 
which is the areal density of a target. One can expect that the applicability of the Moliere's 
formula is restricted by the low density. However, the experimental results obtained for 
small scattering angles show that the deviations from the Moliere's formula are small for 
all data available. In the present paper we explain such surprising behavior calculating the 
leading bulk-density correction to the Moliere's formula. 

We start with the expression for the small-angle scattering amplitude. This expression 
has been obtained in Ref. |^ in the quasiclassical approximation with the first correction 
taken into account. The applicability of this approximation is provided by small scattering 
angles and high energy e of the particle, e m (m is the particle mass, the system of units 
with ^ = c = 1 is used). This amplitude has been obtained for arbitrary localized potential 
without requirement of its spherical symmetry. As known, the quasiclassical wave function 
has a much wider region of applicability than the eikonal wave function. However, as it 
was shown in Ref. [9|, the scattering amplitude obtained with the use of the quasiclassical 
wave function coincides with that obtained in the eikonal approximation, see also Ref. (lo| . 
Using the quasiclassical scattering amplitude with the first correction taken into account, 
we calculate the corresponding cross section and average it over the positions of atoms in 
the target. Dividing this cross section by the area of the target, we arrive at the angular 
distribution The leading term of this distribution coincides with the Moliere's formula. 
The correction depends not only on the areal density riL of the target but also on the bulk 
density n alone. We discuss the magnitude of the correction for different target parameters 
and scattering angles. 
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II. DIFFERENTIAL PROBABILITY 



Let us direct the z axis along the initial momentum p of the particle so that r = -zp/pH-^. 
The small-angle high-energy scattering amplitude in the localized potential V {z, $,) has the 
form 



— I dxxA^V{x,^) 



OO X 



dx dyy {V^V{x,^))-{V^V{y,0) 



-OO —OO 



(1) 



where q = p' — p, p' is the final momentum, /C (^) = J_^dxV (x,^), = d/d^, and 
= V|. The second term in braces in Eq. 1^ corresponds to the correction. For q 7^ 0, 
the unity in the leading term can be omitted. The differential cross section, corresponding 
to the amplitude / and having the same accuracy as Eq. reads 



-iq-(€i-€2)g-i[/c($i)-/c(^2) 



1 + 1 y dxxA^,V{x,^,)-'- J dx J dyy{V^,V{x,^i))-{V^,V{y,^,))y (2) 

—00 —00 —00 

The total potential of atoms in the target has the form 



V (r) = '^u{r - Ti) 



(3) 



where u (r) is the potential of individual atom, which we assume to be spherically symmetric. 
We perform the averaging over the atomic positions using the prescription 



(/)=/n^/. 



(1) 



corresponding to the dilute gas approximation. As a result, we obtain 



3 h 



(5) 
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where 



^1 = y ^exp{-i[x(pi) -x(pi -p)]}, 

P^ = \j ^exp{-z[x(pi)-x(pi-p)]}[(Vp,x(Pi))-(Vp,x(Pi))+^Ap,x(Pi)] 

OO X 

+ j dx j dyy j ^exp{-i [x (pi) - x (pi - p)]} (Vp^n (a;, pi)) • (Vp^n (y, pi)) 
dpi dp2 



— OO — OO 

OO OO 



X 



dx J dy 

OO — OO 

L _ {x - yY 
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S S 



exp {-i [x (pi) - X (pi - p) + X (P2) - X (P2 - P)]} 



4L 



(VpiM(x,pi)) ■ (Vp2M(?/, P2)) . 



(6) 



Here x (p) = c^^^m (x, p), so that /C (p) = x{p ~ Pi)- Then we pass to the limit 
N, S ^ 00 and N/S = riL = const. In this limit, 



= (1 + / ^[exp{-^[x(Pi)-x(Pi-p)]}-l]) 
-nL j dpi (1 - e-*[^(''i)-^(''i-'')]) 



AT 



exp 



Substituting Eqs. (I6|) and ([7]) into Eq. ((51), we finally obtain 



(7) 



dW 

In 



dpe exp <^ -nL / c/pi [1 - cos (x (pi) - X (Pi - p))] 



X 1 1 - ^ y rfx y c/pi sin (x (pi) - X (Pi - P)) Pi ■ VpjM^ (x, pi) 



dpi cos (x (pi) - X (Pi - P)) PiX (Pi 



X j dp2 sin (x (P2) - X (P2 - p)) Vp,x (P2) j- (8) 

At the derivation of this formula we have used the integration by parts. As a result, all 
terms proportional to L in F2 and F3, Eq. (JG]), vanished. It is convenient to rewrite Eq. 
(E]) in another form. The differential in the momentum transfer Q cross section of high- 
energy scattering off one atom, calculated in the quasiclassical approximation with the first 
correction taken into account, satisfies the relation 

I d'Q (1 - e^^-^) ^ = Jdpi[l- cos (x (Pi) - X (Pi - P)) 

00 

+ ^ J rfxsin(x(pi) -x(pi -P))pi ■ VpiM^(x, pi) . (9) 



Using this relation, we obtain with the same accuracy as Eq. (I8|) the following form 



dW 



(27r) 



X 1 



2 / dpe exp 



da 
dQ 



dpi cos (x (pi) - X (Pi - P)) PiX (Pi 



X j dp2 sin (x (p2) - X (p2 - p)) Vp^x (P2) |. 



(10) 



The leading term in Eq. (|TOl) . corresponding to unity in braces, coincides with the 



Moliere's formula. The correction describes the effect of the bulk density of the target 
and has not been known so far. 



III. DISCUSSION 

Let us discuss the magnitude and the structure of the correction obtained. At fixed areal 
density nL (the number of target atoms per unit area), the correction behaves as n (or 
L~^), and increases when L decreases. Estimations show that the relative magnitude of the 
correction is the largest when the main contribution to the integral over p in Eq. (fTOll comes 
from the region p ~ a, where a is the screening radius of atom, a ~ a^Z^^/^, is the Bohr 
radius, Z is the nuclear charge number. This condition is fulfilled when q ~ nLa, where q 
is the momentum transfer. In this case, the correction has the relative order 

where a = 1/137 is the fine-structure constant. Using the estimates 

{eay^ < {may^ ~ aZ^^^ < 1, 
Zana^ < Zaa^^ (a^Z^^/^)^ = « < 1, 

we obtain 

e 

The upper bound for 5 grows with R. However, when i? ^ 1, both the leading term and the 
correction are suppressed by the factor exp [—hR] , where 6 ~ 1 is some numerical constant. 
Therefore, in the whole region interesting from the experimental point of view R is not too 
big so that 5 is very small. 



5 



To conclude, we have calculated the volume density correction to the Moliere's formula 
and estimated its magnitude. This correction turns out to be very small for all reasonable 
values of parameters. Therefore, the Moliere's formula remains very accurate even for high 
density of the target. 
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